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Abstract 

We develop a mathematical model for the energy landscape of poly¬ 
hedral supramolecnlar cages recently synthesized by self-assembly [H]- 
Our model includes two essential features of the experiment: (i) geom¬ 
etry of the organic ligands and metallic ions; and (ii) combinatorics. 

The molecular geometry is used to introduce an energy that favors 
square-planar vertices (modeling Pd^’*’ ions) and bent edges with one 
of two preferred opening angles (modeling boomerang-shaped ligands 
of two types). The combinatorics of the model involve 2-colorings of 
edges of polyhedra with 4-valent vertices. The set of such 2-colorings, 
quotiented by the octahedral symmetry group, has a natural graph 
structure, and is called the combinatorial configuration space. The 
energy landscape of our model is the energy of each state in the com¬ 
binatorial configuration space. 

The challenge in the computation of the energy landscape is a com¬ 
binatorial explosion in the number of 2-colorings of edges. We describe 
sampling methods based on the symmetries of the configurations and 
connectivity of the configuration graph. When the two preferred open¬ 
ing angles encompass the geometrically ideal angle, the energy land¬ 
scape exhibits a very low-energy minimum for the most symmetric 
configuration at equal mixing of the two angles, even when the average 
opening angle does not match the ideal angle. 
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1 Introduction 

The central goal of the field of synthetic self-assembly is to create ordered 
assemblies of building blocks by mimicking the essentials of biophysical pro¬ 
cesses of self-assembly. In supramolecular chemistry, the building blocks are 
molecules, and a basic problem is to synthesize supramolecular cages, to 
shed light on the self-assembly of viral capsids 13 El [I]. Such cages may 
be used to contain or catalyze chemical reactions with great precision. For 
example, molecular flasks that stabilize highly reactive compounds, such as 
phosphorus, have been synthesized by self-assembly [6]. While the chemistry 
of these examples differ, all of them share fundamental geometrical features. 
The cages are Platonic or Archimedean polyhedra, and their molecular de¬ 
sign relies on the decomposition of the polyhedra into simpler geometric 
units, and a careful choice of molecules whose reactive units can be ide¬ 
alized as such units. For example, in the experiments reported in [^, the 
fundamental units are hexagonal ‘molecular tiles’ that assemble into (quasi)- 
truncated octahedra, each molecular tile forming one hexagonal face of the 
truncated octahedron. 

Our primary interest in this article is to develop a mathematical model 
that sheds light on the self-assembly of organometallic molecular cages syn¬ 
thesized in Fujita’s lab [9]- These cages may be idealized as Platonic or 
Archimedean solids with 4-valent vertices. They are realized in experiments 
by organometallic complexes consisting of 4-valent palladium (Pd^^) ions 
at vertices linked by ‘boomerang-shaped’ organic ligands. All such molec¬ 
ular cages have the general chemical formula M„L 2 n- The subscript of M 
indicates the number of palladium ions at the vertices, and the subscript of 
L denotes the number of organic ligand molecules making up the edges of 
the polyhedron. The symmetry of the vertices forces n = 6, 12, 24, 30 and 
60. Several of these ligand polyhedra were observed in experiment, including 
the octahedron (n = 6), cuboctahedron (n = 12), and rhombicuboctahedron 
(n = 24) [9]. 

A particularly interesting effect was observed when two different ligand 
molecules with differing bend angle were mixed. The concentration of the 
two different organic ligand molecules, types A and B say, served as a con¬ 
trol parameter. Separately, these ligands formed cuboctahedra and rhom- 
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bicuboctahedra of the form Mi 2 d .24 and M 24 -B 48 respectively. In a mixture, 
a variety of supermolecules with the general chemical formula MnA 2 n-mBm, 
0 < m < 2n, and n = 12 or n = 24 are theoretically possible. Further, for 
1 < m < 2 n — 1 these polyhedral supermolecules typically have many dis¬ 
tinct isomers, corresponding to different geometric arrangements of the A 
and B ligands. In experiment, as the concentration ratio was varied, a sharp 
transition was observed from a solution in which only cuboctahedra formed 
to one in which only rhombicuboctahedra formed. More precisely, while a 
solution consisting of all isomers MnA 2 n-mBm with n = 12 or n = 24 and 
0 < m < 2 re is possible in theory, in each experiment with a hxed con¬ 
centration ratio of A and B ligands, the equilibrium solution consisted of 
only cuboctahedral isomers (n = 12 ) or only rhombicuboctahedral isomers 
(n = 24). The sharpness of this phase transition was described as emergent 
behavior by Fujita and co-authors [9]. In order to understand what drives 
this phase transition, we develop a ‘minimal’ model that combines the ge¬ 
ometry of isomers with a phenomenological energy for deviations from the 
ideal bend angle. We then compute an energy landscape for isomers with 
the general chemical formula MnA 2 n-mBm, n = 6, 12 and 24, 0 < m < 2n. 

Despite the apparent specificity of our model, the ideas presented here 
may be naturally adapted to the other experiments. More broadly, the 
development of a theoretical understanding of self-assembly is a deep fun¬ 
damental exercise, rich in biological, mathematical and physical ideas. Our 
work is an instance of small, but growing, mathematical literature on syn¬ 
thetic self-assembly that emphasizes ideas from discrete geometry and sta¬ 
tistical physics ia0ii!. These models typically include three aspects. The 
first is the identihcation of a conhguration space that idealizes the set of 
intermediate structures that lie between the building blocks and the assem¬ 
bled product. The second is to understand the free energy landscape; in 
particular, to understand which of the states in the configuration space are 
energetically favorable. The third is a description of the kinetics of self- 
assembly. In this article, we focus on the first two aspects of this procedure, 
since (as will become clear) determining the energy landscape requires in¬ 
tensive computation, and is of independent interest. 

The rest of this article is organized as follows. We describe the configura¬ 
tion space and energy in the next section. This is followed by a description 
of the combinatorial explosion, and a symmetry-based sampling scheme. 
Finally, we describe our numerical results and conclusions. 
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2 The configuration space and energy function 


We are interested in molecules of the form MnA 2 n-mBm, with n vertices 
corresponding to the ions M, 2n — m ligands of type A, and m ligands 
of type B. Our model involves: (i) simple geometric combinatorics, (ii) 
a phenomenological energy function, and (iii) sampling schemes for fast 
computation, (i) For any m, the number of distinct isomers is the number 
of distinct 2-colorings of the edges of the polyhedra modulo the action of 
the symmetry group of the polyhedron, (ii) Our energy function penalizes 
the geometric distortion of the polyhedra caused by ligand molecules with 
differing bend angles, (iii) A calculation of the full energy landscape is 
impossible for n = 12 and n = 24 because of a combinatorial explosion. We 
use symmetry and heuristics based on n = 6 to sample the energy landscape. 

2.1 Colorings and configurations 

More formally, let "P be a polyhedron whose faces, edges and vertices are 
denoted We assume that P is one of the following: the octahe¬ 

dron; cuboctahedron; or rhombicuboctahedron. All these polyhedra, shown 
in Fig. have regular 4-valent vertices and their symmetry group is the oc¬ 
tahedral group, 0/i [2]. A 2-coloring of edges, or coloring for short, is a map 
c : £ —7- {0,1}. The set of colorings is denoted C. It is convenient to view 
the colorings ‘physically’: we identify the polyhedron P with its standard 
embedding in space, and paint the edges of P blue (0) or red (1), for ligands 
of type A and B respectively. We say that two colorings ci and C 2 are equiv¬ 
alent under Oh, written ci ~ C 2 , if there exists g € Oh such that gci = C 2 . 
We define a configuration to be an equivalence class of colorings under the 
relation ~. The quotient set CjOh is called the combinatorial configuration 
space, or simply configuration space. It is the set of all distinct colorings 
modulo the symmetry of the polyhedron. 

The enumeration of CjOh is a classic problem in discrete group theory 
(we follow 0). Let OhC denote the orbit of a coloring c. We say that 
the degeneracy of a coloring c is the size of the orbit \Ohc\. Clearly, the 
degeneracy is the same for all colorings in the equivalence class [c]. Thus, 
it is meaningful to speak of the degeneracy of [c]. Let Oc C Oh denote the 
stabilizer subgroup of c £ C. We call the size of this subgroup, |oc|, the 
symmetry number of c, and denote it by s. By the orbit-stabilizer theorem 
and Lagrange’s theorem, for every coloring c £ C, 

\Ohc\ X \o,\ = \Oh\ = 48. (1) 
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Figure 1: 4-valent polyhedra with Oh symmetry, a) octahedron; b) 
cuboctahedron; c) rhombicuboctahedron. (d-f) Ideal supramolecular 
cage embeddings. Here, the Pd^"*" ions are represented by large black 
vertices, while the ‘boomerang’-shaped organic ligands are represented by 
small gray vertices (‘elbows’), each with two arms; in (d), the straight edges 
of the octahedron are included only for comparison. Note that the geometry 
at the ions is square-planar, (g-i) Frustration. Minimum-energy embed¬ 
dings when the preferred bend angle of the edges is 9 q < 9^. g) 9o = 81°, 
(9* = 90°; h) 9o ~ 105.5°, « 117.2°; i) 6»o « 120.7°, 0* ss 134.1°. Bright¬ 

ness of color indicates the energy contribution from each vertex and each 
‘elbow’. Note that the geometry at the ions is deformed from square-planar. 
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Two colorings ci and C 2 are equivalent if and only if they lie on the same 
group orbit. Thus, the size of the configuration space, \C/Oh\, is simply 
the number of distinct orbits. Let denote the set of colorings fixed by 
g G Oh, i.e. C® = {c G C \gc = c}. By Burnside’s lemma. 


\C/Oh\ 


1 


E i^^i- 

g&Oh 


( 2 ) 


While this calculation of the size \C/Oh\ is well-known, our interest lies 
in an explicit description of each equivalence class [c] G CjOh- We compute 
CjOh explicitly for the octahedron and cuboctahedron as described below. 
We are unable to explicitly compute CjOh for the rhombicuboctahedron 
because of a combinatorial explosion. Note that \C/Oh\ is bounded below 
by 




2l^l 

4 ^' 


(3) 


The octahedron has 144 unique configurations (the lower bound is 85) and 
the cuboctahedron has 352,744 unique configurations (the lower bound is 
349,525). For the rhombicuboctahedron, the same bound yields 2^®/48 ~ 
5,864,062,014,806, which is too large to enumerate explicitly. Thus, one 
of the main computational challenges in our work is to obtain a realistic 
understanding of the energy landscape, despite the combinatorial explosion 

oi\C/Oh\. 

Finally, note that the conhguration space CjOh inherits a natural graph 
structure from C. We say that two colorings Cm and Cm+i with m and m -|-1 
‘red’ edges are neighbors in C if they differ at a single edge. If Cm and Cm+i 
are neighbors, it is clear that gcm and gcm+i are also neighbors for each 
g ^ Oh- Thus, every coloring in the orbit OhCm is the neighbor of at least 
one coloring in the orbit OhCm+i and it is natural to say that [cm] and [cm+i] 
are neighbors in CjOh- 


2.2 A phenomenological energy for molecular distortion 

Each configuration [c] G CfOh determines the unique combinatorial struc¬ 
ture of a supramolecular conformation. An energy landscape is a map 

F : C/Oh ^ M. 

We define a phenomenological energy by penalizing the distortion both 
of the edges and of the vertices from their preferred geometry. The edges 
are modeled as two straight arms, joined at an angle at the ‘elbow’ of the 
boomerang (Fig. [^). Let X G m3|£^|-i-3|V| (denote the coordinates of the 
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Figure 2: Geometry of a) edges, comprising an ‘elbow’ and two arms; and 
b) 4-valent vertices. 


elbows of the ligands and of the vertices of the embedded polyhedron. We 
introduce energy terms which account separately for the distortion from 
a preferred length for each arm of an edge {Earms)] distortion from the 
preferred bend angle at each elbow {Ef,ibow)\ and distortion from square- 
planar geometry at each vertex {E^ertex)'- 

G{X- [c]) = {Earms{t{X)) + E^MX))) + Y, E^ertexiH^)) ■ (4) 

ce£ DSV 

Here the notation e(X) and t)(X) implies that the energy depends on the co¬ 
ordinates of each edge and vertex respectively. Finally, the energy landscape 
F([c]) is defined by 

F([c]) = infG(X;[c]), (5) 

Ji. 

where the infimum is taken over all embeddings of V into In our nu¬ 
merical experiments, the energy G is minimized through a zero-temperature 
Monte-Carlo annealing. This annealing does not reach the exact energy min¬ 
imum, but gives an upper bound on the minimum-energy which we believe 
to be typically within a few percent of the actual infimum. This fact may 
be verified for the minimum energy configurations shown in Fig. EKd-f). 

We now describe the terms in the energy in local coordinates in more 
detail. We introduce an energy cost to deviations of the length of the arms 
from one: 

Earm = Ea{{h - l.O)^ + {h - l.O)^), (6) 

where li and I 2 are the lengths of the two arms. We choose Ea ^ Ee,Ep, E^ 
(defined below) to make the arms stiff; this has the effect that in general 
in the minimum-energy embedding, the arms are of length one and do not 
contribute signficantly to the total energy. 
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We also introduce an energy cost to deviations of the angle of the elbow 
from its preferred opening angle: 


Eelbow = Ee{6 - 0o)^ (7) 

with 0 the actual opening angle and 6q the preferred opening angle (which 
depends on the type A or B). Thus this term is the only term in the 
energy which depends on the configuration as well as on the embedding. 
The coefficient sets the scale of the energy contribution of the edges. 

At each vertex, we introduce an energy cost for deviations from square- 
planar geometry (see Fig. [^): 


E.. 


vertex 


— E^, 


- 7r/2f + EpY,i^i - ^/ 2 )' 
(i,j) * 


( 8 ) 


with the first sum over pairs of neighboring edges i and j and 9ij the angle 
formed by the ‘arms’ of those edges; and the second sum over edges i and 
the angle formed between the arm of that edge and ‘plane perpendicu¬ 
lar’, defined as the average of the cross-products of the pairs of neighboring 
edges. The first sum pushes the edges to all be perpendicular to one an¬ 
other (square), while the second sum pushes the edges to be co-planar. The 
overall coefficient sets the scale of the energy contribution of the vertices 
relative to the edges, while the coefficient Ep sets the scale of the energy 
contribution of the ‘planar’ aspect of the vertices relative to the ‘square’ 
aspect. 

For each polyhedron, there is an ideal angle 0* such that there is a zero- 
energy state when all edges have Oq = 6^. We calculate these ideal angles 
analytically for each of the three polyhedra considered. Fixing one vertex at 
the origin, we write explicitly the positions A, B, C, D of its four neighbors 
in a regular polyhedron (ordered cyclically so that B and D are each closer 
to A than is C). We then solve for E, F, G, H, the positions of the elbows, 
which satisfy the following conditions: 


|E| = |F| = |G| = |H| 

G = -E , H = -F 
EF = EH = FG = GH = 0 

A - E| = |E| , |B - F| = |F| , |C - G| = |G| , |D - H| = |H 
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These relations ensure that the geometry at the vertex is square-planar, and 
that every arm of an edge is of the same length. Examining the triangle 
OEA we see that the opening angle 9 is then given by 

(' ■ iw) 

Thus if the preferred angle of every edge is set to 0*, the embedding with this 
geometry at every vertex of the Platonic or Archimedean polyhedron will 
give a zero-energy configuration; these embeddings are shown in Fig. [^d)- 
(f). The values derived for the ideal angles of the three polyhedra are: 


,octahedron 
jCuboctahedron 


,rhombicuboctahedron 


= 90°, 
= cos' 


= cos 


(^] 

(25- 20V2 

\33 - 20 V 2 


PS 117.2°, 

Ri 134.1°. 


Note that these values for the cuboctahedron and rhombicuboctahedron 
correct the approximate values of 120° and 135°, respectively, given in [9]. 

However, in experiment this value can only be chosen approximately, 
and the angle Oq is different from the ideal value. This causes frustration, as 
shown in Fig. [^g)-(i), and the minimum-energy state balances competing 
energies. In these examples, the total energy of the vertices is between 1.8- 
2.0 times the total energy of the edges. All of the energy contributions are 
necessary to observe frustration. This can be seen by counting the degrees 
of freedom. For example, if the planarity of the vertices is relaxed, the edge 
‘elbows’ can rearrange themselves so that they satisfy their preferred angles 
at the expense of the square-planar geometry at the vertices. Only when all 
four energy terms are in competition do we have an interesting system. It 
is in this sense that we consider our energy model to be minimal. 

Finally, we define the average preferred angle for mixtures of two ligands. 
The preferred angle used in [9] is: 

9q{ca, cb) =-,- {ca6q{A) + cb9q{B)) , (10) 

CA + CB 

where ca and cb denote the concentrations of ligands A and B. We will 
also use 9q with no arguments to mean 

eo^eoic,c) = ^{9o{A) + eom (ii) 
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and 0o{m) in an isomer MnA 2 n-mBm to mean 

9o{m) = (9o(2n — m, m) = — ((2n — m)9o(A) + m9o{B)) (12) 

2n 

3 Symmetries and Sampling Strategies 

As discussed above, our goals are twofold: 1) to enumerate the configura¬ 
tion space Cl Oh] and 2) to calculate the energy landscape F{C/Oh)- The 
computational difficulty of these tasks depends strongly on the size of the 
polyhedron V. When V is the octahedron, we are able to fully enumer¬ 
ate the configuration space and its graph structure, and compute the entire 
energy landscape. When V is the cuboctahedron, we fully enumerate the 
conhguration space and its graph structure, but we sample the energies of 
only a fraction of the conhguration space to compute a partial energy land¬ 
scape. When V is the rhombicuboctahedron, even the conhguration space 
is intractable. We enumerate only those conhgurations c with non-trivial 
symmetry number s, and sample the energy of only a fraction of these con¬ 
hgurations. In what follows, we focus mainly on the sampling strategies 
used to obtain partial, but informative, descriptions of the conhguration 
space and energy landscape for the two larger polyhedra. 

The octahedron is the most straightforward. With only 2^^ = 4096 
colorings to consider, we construct every coloring, and calculate the equiv¬ 
alences under Oh to obtain a complete list of conhgurations. There are 144 
conhgurations, of which 94 have non-trivial symmetry groups (see Table [^. 
We also calculate explicitly the graph structure, identifying the neighbor 
relations between conhgurations. We then calculate the energy of every one 
of the 144 conhgurations by solving the minimization problem ([^. 

We build the graph Of Oh for the cuboctahedron by induction on the 
number of ‘red’ edges, denoted m. There is a single conhguration which 
has m = 0, that is, all edges are ‘blue’; this provides the base case. As 
our induction step, we use the subset {CfOh)m of conhgurations with m 
‘red’ edges to construct the set {C/Oh)m+i as follows. For each element 
[c] G {CjOh)m, we consider a representative coloring c G [c]. This coloring 
has neighbors {c'} C Cm+i, obtained from c by changing the color of a single 
edge from ‘blue’ to ‘red’. We compute the orbit Ohc' for each neighbor c', 
and compare each element of Ohd against the conhgurations already seen in 
{C/Oh)m+i- If the conhguration has already been seen in {C/Oh)m+i-, we 
make note that [c] and [d] are neighbors. If the conhguration has not been 
seen before, we update our list of representatives in {C/Oh)m+i to include 
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symmetry 

number 

octahedral 

configurations 

cuboctahedral 

configurations 

rhombicuboctahedral 

configurations 

48 

2 

2 

4 

24 

0 

1 

4 

16 

2 

2 

28 

12 

2 

6 

29 

8 

6 

23 

644 

6 

4 

35 

716 

4 

21 

211 

46,991 

3 

1 

44 

16,015 

2 

56 

5,943 

>40,648,786 

1 

50 

346,477 

PS 6 X 10^2 

total 

144 

352,744 

PS 6 X 10^^ 


Table 1: Enumeration of configurations by symmetry number for each poly¬ 
hedron. 


[d] and make note of the neighbor relationship. This process ends when we 
find the single configuration with m = 2n = 24. Thus we do not enumerate 
every possible coloring, but we nonetheless identify all configurations and 
compute the connectivity of the graph CjOh- We find a total of 352,744 
configurations (see Table [^. 

The computation of the energy for each configuration [c] involves a nu¬ 
merical solution of the minimization problem ([^. Rather than pursue this 
computation for each of the 352,744 configurations, we sample a fraction of 
the configurations and compute a partial energy landscape. We use a sam¬ 
pling strategy which is designed to identify the full range of energies at each 
value of m, that is, to find the lowest and highest energy configurations at 
each m. The strategy exploits two heuristic observations about the energy 
landscape of the octahedron: 

1. For a given m, the highest and lowest energy configurations tend to 
have nontrivial symmetry number s; 

2. High energy configurations tend to have high energy neighbors in the 
graph. Likewise, low energy configurations have low energy neighbors. 

Our sampling strategy consists of two steps. First, we calculate the 
energy for all configurations with ‘extreme’ values of m, here defined as 
m < 4 or m > 20 (630 configurations total), and for all configurations 
with symmetry number s > 3 (324 configurations total, although there is 
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overlap with the 630 extreme-m configurations). Finally, we calculate the 
energies of about 2 % (but at least 100 ) of the remaining configurations 
for each value of m such that 5 < m < 19, chosen as follows. About ^ 
to ^ of these configurations are identified as the neighbors of the highest 
and lowest energy configurations known at m — 1 and m + 1. The other 
roughly half of these configurations are sampled uniformly at random from 
the remaining configurations. In order to maximize the information we have 
about neighbors of high and low energy configurations, we carry out this 
calculation sequentially for different values of m, starting with m = 5 and 
m = 19 and working our way to less extreme values, ending with m = n = 
12 . 

For the rhombicuboctahedron, the combinatorial explosion in the num¬ 
ber of configurations makes it intractable to even enumerate the approxi¬ 
mately six trillion configurations. We focus instead on configurations with 
non-trivial symmetries, motivated by the heuristics above. 

A configuration with a non-trival symmetry number s > 2 is fixed by 
at least one non-identity element g G Oh. We generate such ‘symmetric’ 
configurations by first enumerating the set of colorings fixed by 5 , , for 

each g ^ Oh, g ^ e as follows. We consider the orbits of the edges of the 
polyhedron under the subgroup Sg = (g) C Oh generated by g. In order for 
a coloring c to be fixed by g, i.e. gc = c, the edges in each orbit SgC, z G S 
must be of the same color. Thus, is constructed as the set of 2-colorings 
of edge orbits, SgZ. We classify the forty-seven non-trivial elements of Oh 
into the following five subsets: 

(i) There are 20 elements of Oh with order > 4. Each of these elements 
has a power g°' which is of order 2 or 3. Since any coloring fixed by g 
is also fixed by g°‘, C , and we need only enumerate for g of 
order 2 or 3 to get all configurations. 

(ii) There are 8 elements of Oh with order 3; while these fix different 
colorings, they generate the same equivalence classes in Of Oh- Thus 
we need enumerate for only one of these 8 elements. 

(iii) There are 10 elements of Oh with order 2 which have the property that 

= 2 for all e G T. Of these only three produce different config¬ 
urations, while the remaining seven fix colorings which are equivalent 
to colorings already identified. 

(iv) There are 6 elements of Oh with order 2 which have four orbits for 
which |S'ge| = 1 , that is, these elements map four edges onto them- 
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selves. As with the elements of order 3, these produce the same con¬ 
figurations, and so we enumerate for only one. 

(v) There are 3 elements of Oh with order 2 which have eight orbits for 
which |Sge| = 1, that is, these elements map eight edges onto them¬ 
selves. These elements represent reflections across the plane bisecting 
a ring of square faces. Again, we enumerate for only one of them. 

Having enumerated for several elements g, we obtain a list of colorings 
which contains at least one representative of [c] for each [c] G C /Oh such that 
s > 2. We then compare these colorings and their equivalences to enumerate 
all [c] with s > 3 (64,431 configurations total; see Table[^. There are tens of 
millions of configurations with symmetry s = 2; we have not enumerated all 
of these, but have enumerated those with extreme values of m {m < 8 and 
fn > 40), along with those at values of m relatively prime to 48 which have 
no higher-symmetry configurations (m = 11,13,17,19, 23, 25, 29, 31,35, 37). 

We then calculate energies for configurations sampled as follows. We 
calculate the energy for all configurations with symmetry s > 6 (1425 con¬ 
figurations total); about 2.5% (but at least 40 at each m) randomly selected 
configurations with s = 4; about 2.5% (but at least 40) randomly selected 
configurations with s = 3; and for m < 8, m > 40, or m relatively prime to 
48, about 0.02% (but at least 200) configurations with s = 20 This gives 
us a total of roughly 10,000 energies calculated for each of the plots in the 
next section - a small fraction of the 6 trillion configurations, but a sam¬ 
pling which we believe captures the most interesting aspects of the energy 
landscape. 


4 Results 


In most of the examples below, we set the energy parameters Ea = 10.0; 
Ep = 1.0; and Eg = Ey with 1.0 < E^. < 3.8, with the precise value typi¬ 
cally chosen to normalize the energy of the m = 0 configuration to 1.0. In 
Section 4.5 we describe results with E^, / Ey. 


4.1 Energy landscape when 9 q = 9 ^ 

In Fig. 1^ we present sample energy landscapes for the octahedron, the 
cuboctahedron, and the rhombicuboctahedron, when the average preferred 

^In Fig. [^, we calculate the energy for configurations with s = 3 only for m < 8 and 
m > 40, and for a larger proportion of configurations with s = 2 for m < 8, m > 40, 
m= 11,13,17,31,35,37. 
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angle is the same as the ideal angle, that is, 0 * = 0 o = 5 (^o(^) + ^o(-B)). 
(This choice of parameters was used as a heuristic design principle in 
as discussed below, the broad features of this energy landscape are largely 
insensitive to these choices of Oq{A) and 6q{B).) We observe that when 
n and m are both held fixed, there is nonetheless a strong variation in 
the energy of all isomers with chemical formula MnA 2 n-mBm- The overall 
maximum energy is found at m = n, as is the minimum energy for both the 
cuboctahedron and the rhombicuboctahedron; this minimum energy appears 
to be zero, within the error of our minimization. For the octahedron, the 
minimum energy is 1.0 at m = 0 or m = 2n, but there is a low-energy 
configuration at m = n with only slightly higher energy. 

4.2 Extremal energy configurations have high symmetry 

This low-energy configuration at m = n is also the isomer with the largest 
symmetry number. While this may seem unsurprising, it must be recognized 
that there are also high-energy configurations with high symmetry number! 

More generally, in all of the energy landscapes shown in Fig. for any 
given m, the configurations with non-trivial symmetry number tend to span 
or nearly span the range of energies. For the octahedron, for which the 
energy landscape is completely calculated, the extremal energy configura¬ 
tions at each m have symmetry number s > 2; indeed, for m 7 ^ 5, 7, the 
lowest-energy configuration has s > 4. For the cuboctahedron, we have cal¬ 
culated only a partial energy landscape, with our sampling biased toward 
high-symmetry configurations, but nonetheless we find similar results. 

Some insight into why symmetry plays such an important role may be 
gained by considering the structure of these extremal-energy configurations, 
shown in Fig. for m = n for each polyhedron. Consider first the low- 
energy isomers: in all three polyhedra, at every vertex meet two edges of 
type A and two edges of type B in a “'cis''' arrangement such that the 
edges of type A are adjacent to one another, and opposite the edges of type 
B. This arrangement allows the frustration to be resolved by canting the 
plane of the ion (the resolution is not complete for the octahedron, since 
canting fully would require stretching the stiff ‘arms’). Heuristically, then, 
for other values of m, we expect that configurations with many vertices with 
this energetically favorable configurational structure will have low energies, 
and the similarity of vertices tends to yield high symmetry number. In 
the case of the high-energy isomers, those for the cuboctahedron and the 
rhombicuboctahedron show very similar configurational structure: they are 
composed of 6 - or 8 -element closed rings which are either entirely of type 
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Figure 3: Energy landscapes and partial energy landscapes for 

a) the octahedron; b) the cuboctahedron; and c) the rhombicubocta- 
hedron. For each polyhedron, 6q{A) and 9q{B) are chosen so that 
0o(^) <e^< eo{B), and Oq = 9^: a) 9o{A) = 81°, 9o{B) = 99° (0* = 90°); 

b) 9o{A) « 105.5°, 9o{B) « 128.9° (0* « 117.2°); c) 0o(7l) ~ 120.7°, 

9q{B) ~ 147.6° (0* ~ 134.1°). The energies are normalized so that the con¬ 
figuration with m = 0 has energy 1.0. For (a), the entire energy landscape 
is shown; for the partial landscapes in (b) and (c), the choice of sampled 
energies is described in the text. Color indicates the symmetry number of 
the configuration. 
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A or entirely of type B, so that at every vertex meet either four edges of 
the same type, or two edges of each type in a “trans” arrangement such 
that edges of the same type lie opposite one another. The largest energy 
contributions are from vertices at which four edges of type B meet, followed 
by those at which four vertices of type A meet. Again, configurations with 
many vertices with these energetically unfavorable structures will have high 
energies. 

4.3 Low and high energy configurations are neighbors in the 
graph structure 

For the octahedron, we obtain a complete description of the configuration 
space as a graph, as shown in Fig.j^ We observe that low (respectively high) 
energy configurations are usually neighbors of other low (resp. high) en¬ 
ergy configurations. This feature of the energy landscape may be explained 
heuristically as follows. Starting from a low (resp. high) energy config¬ 
uration with several favorable (resp. unfavorable) vertex arrangements as 
described above, changing the type of a random edge will leave many of these 
vertex arrangements unchanged, resulting in a modest change in energy and 
giving another low (resp. high) energy configuration. 

4.4 Energy landscape when Oq ^ 9^ 

We find that the shape of the energy landscape is robust for variations of 
do{A) and do{B). When 6o{A) < 9^ < Bo{B), but 0* / 0o = U^o{A) + Bo{B)), 
the energy spectrum becomes tilted as shown in Fig. [^, with 

E{m = 2n) ^ \Bo{B) — 

E{m = 0) \Bo{A)-B,\‘^' 

Importantly, the lowest energy is still found at even mixing, m = n, even 
though Bo{m = n) ^ B^; the energy is nonzero, but still small compared to 
the energy at m = 0 and m = 2n. We conclude that even if Bq ^ B^, the 
energetically favorable arrangements of type A and type B edges described 
above allow for much of the frustration to be relieved. 

The lowest energy configuration instead occurs for m = 0 when B^ < 6o{A) < Bo{B), 
and the energy increases roughly linearly with m, as shown for the cuboc- 
tahedron in Fig. [6>. Isomers with a fixed m still exhibit a range of energies, 
but the variation with m is similar to the m-dependent energy contribution. 
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Figure 4: Extremal energy configurations for m = n. For clarity, the 
edges are drawn as straight rather than showing the full embedded cage 
structure. The preferred angles are as in Fig. blue edges are of type A 
with 6o{A) < 0*, while red edges are of type B with 6o{B) > 9^. Bright¬ 
ness of the color indicates the energy contribution from each vertex and 
each edge; the color scale is different in each panel. Note the deformation 
from a regular polyhedron; in particular, blue edges are shorter than red 
edges, a-c) Minimum-energy configurations for a) the octahedron; b) the 
cuboctahedron; and c) the rhombicuboctahedron. d-f) Maximum-energy 
configurations. Note the high symmetry numbers: a) s = 12; b) s = 24; c) 
s = 48; d) s = 2; e) s = 6; f) s = 6. 
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Figure 5: Energy landscape of the octahedron shown with the graph struc¬ 
ture of the configuration space. (The range of the y-axis has been changed 
from Fig. for clarity). Note that low energy configurations tend to be 
connected to other low energy configurations, and likewise high energy con¬ 
figurations to high energy configurations. 

4.5 Independence of energy scales 

There is little variation in the shape of the energy spectrum as the param¬ 
eters Eg, and Ep are varied (we fix Ea = 10.0 > Ee, Ey, Ep). Fig. 
shows the energy spectrum for a rhombicuboctahedron with the same pre¬ 
ferred angles as in Fig. [^, such that Oq = 9^, but with Eg = 2Ey (Ep is kept 
at 1.0). This is nearly indistingnishable in overall shape from Fig. 1 ^, with 
Ee = E^; there is a slight increase in the maximal energies reached but it is 
subtle. Similar results were obtained for the cuboctahedron with E^, = ^Ey] 
again the energy landscape has much the same shape, with in this case a 
slight decrease in the maximal energies (data not shown). 

We also observed that the ratio between the energy of a rhombicuboc¬ 
tahedron with only one ligand type, M 24 T 48 and a cuboctahedron with the 
same ligand type and parameters, M 124 . 24 , is nearly constant over a much 
larger range of energy ratios, varying both E^/Ey, and Ep (data not shown). 
Thus the ratios of the energy parameters in this system seem to be relatively 
unimportant: the competition between the edge energies and vertex ener¬ 
gies results in an embedding which balances those energies in such a way 
that the shape of the energy spectrum is roughly independent of the ratios, 
and this shape is almost entirely determined by the preferred angles. The 
overall energy scale, of course, will be important in any system in which 
temperature plays a role. 
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Figure 6: Variation of the energy landscapes with model param¬ 
eters. a-b) Change of preferred angles. Partial energy spectra for a) 
the rhombicuboctahedron and b) the cuboctahedron with preferred angles 
Oq{A) = 127°, 6o{B) = 149°; these angles are chosen to match those of 
the two organic ligands reported in [9]. For the rhombicuboctahedron (a), 
Oq{A) < 0* < 9q{B) still holds, but 9o{m) = 9^ for m ^ 16; note that the 
lowest energy configuration is still at m = n = 24. For the cuboctahedron 
(b), 0* < 9q{A) < 9q{B)-, note that the energy increases rapidly with added 
type-i? edges, c) Change of relative energy scales. Partial energy spec¬ 
trum for the rhombicuboctahedron with the same preferred angles as in 
Fig. but with Eg = 2E^. Note that the shape of the energy landscape 
is almost unchanged. For (a) and (c), the overall energy scale is chosen so 
that the configuration with m = 0 has energy 1.0. 
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Figure 7: Crossover of cuboctahedral and rhombicuboctahe- 
dral landscapes. Partial energy landscapes of the cuboctahedron and 
the rhombicuboctahedron, both with preferred angles 0q{A) = 121.5°, 
^o('^) — 151.8 , note O^^cuhoct. ^ ^ 0*^rhombicuboct. ^ Blue cir¬ 

cles give the landscape of the rhombicuboctahedron, and red squares that 
of the cuboctahedron. The landscape of the cuboctahedron is modified by 
doubling both m and the energy, giving in effect the total m and total 
energy of two identical cuboctahedra; this allows us to preserve the total 
number of ions and of ligands when we compare the rhombicuboctahedron 
with the smaller cuboctahedron. The energies are normalized so that the 
rhombicuboctahedron with m = 0 has energy 1.0. 

4.6 Crossover of landscapes for different polyhedra 

We compare two different polyhedra to one another as shown in Fig. 
For the cuboctahedron, we scale both m (the number of ‘red’ edges) and 
the energy by two in order to compare energy per ligand. The preferred 
angles are 0o{A) ss 121.5° and 0q{B) k, 151.8°, so that for the cubocta¬ 
hedron, O^^cuboct. < ^o(^) < 0o{B), while for the rhombicuboctahedron, 
Oo{m) = 0*^rhombicuboct. for m K, 20. Note that there is a crossover; for 
m < 4, the cuboctahedron is energetically favored, while for m > 8, the 
rhombicuboctahedron is energetically favored. This is suggestive that as 
the concentrations of type A and type B ligands are varied, a transition 
may be observed from the cuboctahedron to the rhombicuboctahedron, as 
in |9], although in our model, mixing of the two polyhedra is suggested for 
m ~ 6, which was not observed in [^. 


20 







5 Discussion 


We have studied in some detail the energy landscape for a minimal model of 
supramolecular polyhedra formed by mixtures of two types of ligands. We 
have found a strong variation with configuration in the energy of different 
isomers, configurations which have the same number of type-A and of type- 
B ligands. When the preferred angles encompass the ideal angle for the 
polyhedron, that is, 0o{A) < 9^ < 9 q{B), the lowest energies are achieved 
by high-symmetry configurations in which edges meet at cis configurations 
at a vertex, with a very low energy configuration at m = re. While our 
work does not fully explain the emergence seen in [9], for some choices of 
parameters we see a shift in minimum energy states from cuboctahedron 
to rhombicuboctahedron. We also note a striking analogy to another ex¬ 
periment in Fujita’s lab, in which ligands with different lengths but similar 
angles were mixed [inj. For a sufficiently large ratio between the lengths, 
r ~ 2, they observed a mixed yet sorted structure of the form M 12 A 12 B 12 
with a well-defined structure in which ligands of type A also took on a cis 
configuration opposite ligands of type B. Our work shows that similarly 
mixed yet sorted structures are favored when ligands with different angles 
are mixed, although the energy differences are small enough that the sorting 
is not as complete as in [TO] . 

Our model can easily be extended to a system with temperature by 
introducing two additional terms to obtain the free energy. The first term 
is of the form ksT/s, where s is the symmetry number of the configuration; 
at high enough temperature, this term will change the relative free energies 
of configurations at a fixed value of m. The second term is of the form 
ksT^l — where 1 — p and p are the bulk proportions of ligands 

of type A and B, respectively; this term depends only on m and not on the 
configuration, but at high enough temperature will favor configurations with 
compositions near MnA 2 {i-p)nB 2 pn- These temperature-dependent entropic 
terms must be significant in the experimental results of [9], which report 
polyhedra with a range of values of m, indicating that states other than the 
ground state are populated. 

Further improvements to the model would include the introduction of 
chemically relevant values for the energy scales Ea,Ee,E^ and Ep, which 
would provide a more quantitative comparison to the experimental results. 
Also desirable would be a description of the kinetics; we have so far consid¬ 
ered only the energy landscape, but the question of how the edges of different 
types dynamically arrange themselves is crucial to establishing whether the 
system can reach its equilibrium distribution and on what timescale, or 
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whether non-equilibrium states become trapped. 

These results provide insight into the importance of configuration and 
symmetry in the determination of the free energies of self-assembling sys¬ 
tems. Our work also illustrates the importance of exploiting symmetry to 
reduce the size of otherwise combinatorially intractable problems. 
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